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We analyze turbulent velocity signals measured by two 
probes in the atmosphere, both at the height of 35 meters 
but displaced by 40 cm nominally orthogonal to the mean 
wind. Choosing a suitable coordinate system with respect to 
that of the mean wind, we derive theoretical forms for sec- 
ond order structure functions, and fit them to experimental 
data. We show that the effect of flow anisotropy is small on 
the longitudinal component but significant on the transverse 
component. The data provide an estimate of a universal ex- 
ponent from among a hierarchy that governs the decay of flow 
anisotropy with the scale-size. 



Experimental studies of turbulent flows are usually 
limited in the sense that one measures the velocity field 
at a single spatial point as a function of time fl, and 
uses Taylor's hypothesis to identify velocity increments 
at different times with those across spatial length scales, 
R. The standard outputs of such measurements are the 
longitudinal two-point differences of the Eulerian velocity 
field and their moments, termed structure functions: 



S n (R) 



[u(r + R, t)-u(r,t)] 



R 



(1) 



where (•} denotes averaging over time. In isotropic ho- 
mogeneous turbulence, these structure functions are ob- 
served to behave as a power- law in R, S n (R) ~ R^ n , with 
apparently universal scaling exponents ( n j|. 

Recent progress in measurements || and in simulations 
|| begins to offer information about the tensorial nature 
of structure functions. Ideally, one would like to measure 
the tensorial nth order structure functions denned as 

S° 1 -" a '(R) = ([u ai {r + R)-u ai (r)] 
x [u a2 (r + R)- u° 2 (r)} . . . [u a " (r + R) - u a " (r)]> , (2) 

where the superscript «i indicates the velocity compo- 
nent in the direction i. Such information should be useful 
in studying the anisotropic effects induced by all practi- 
cal means of forcing. In analyzing experimental data the 
model of "homogeneous, isotropic small-scale" is univer- 
sally adopted, but it is important to examine the rele- 
vance of this model for realistic flows. One of the points 
of this Letter is that keeping the tensorial information 



helps significantly in disentangling different scaling con- 
tributions to structure functions |1. Especially when 
anisotropy might lead to different scaling exponents for 
different tensorial components, a careful study of the var- 
ious contributions is needed. We will show below that at- 
mospheric measurements contain important anisotropic 
contributions to one type of transverse structure func- 
tions. 

In this Letter we analyze measurements in atmospheric 
turbulence at Taylor microscalc Reynolds numbers of 
about 20,000 0. The data were acquired simultaneously 
from two single-wire probes separated by 40 cm nomi- 
nally orthogonal to the mean wind direction. The two 
probes were mounted at a height of about 35 m above 
the ground on a meteorological tower at the Brookhaven 
National Laboratory. The hot-wires, about 0.7 mm in 
length and 6 /im in diameter, were calibrated just prior 
to mounting them on the tower (and checked immedi- 
ately after dismounting), and operated on DISA 55M01 
constant-temperature anemometers. The frequency re- 
sponse of the hot-wires was typically good up to 20 kHz. 
The voltages from the anemometers were suitably low- 
pass filtered and digitized. The voltages were constantly 
monitored on an oscilloscope to ensure that they did not 
exceed the digitizer limits. Also monitored on-line were 
spectra from an HP 3561A Dynamic Signal Analyzer. 
The wind speed and direction were independently moni- 
tored by a vane anemometer mounted a few meters away 
from the tower. The real-time duration of data record 
was limited by the degree of constancy of the speed and 
direction of the wind. The Kolmogorov scale was about 
0.45. Table 1 lists a few relevant facts the data records 
listed here (this being part of a much larger set). The 
various symbols have the following meanings: U = lo- 
cal mean velocity, v! = root-mean-square velocity, (e) 
= energy dissipation obtained by the assumption of lo- 
cal isotropy and Taylor's hypothesis, rj and A are the 
Kolmogorov and Taylor length scales, respectively, the 
micro-scale Reynolds number R\ = u'X/v, and f s is the 
sampling frequency. 

As a first step we tested whether the separation be- 
tween the two probes was indeed orthogonal to the mean 
wind. To this end we computed the cross-correlation 
function (ui(t + r)u2(t)) where subscripts 1 and 2 re- 
fer to the two probes. If the separation were precisely 
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TABLE I. Basic information about the data analyzed in 
this paper. 



orthogonal to the mean wind, this quantity should be 
maximum for r = 0. Instead, we found the maximum 
at 0.03 sec, implying that the separation was not pre- 
cisely orthogonal to the mean wind. To correct for this 
effect, the data from the second probe were time-shifted 
by 0.03 seconds. This amounts to a change in the ac- 
tual value of the orthogonal distance, and we computed 
this effective distance to be A « 31 cm. Next we tested 
the isotropy of the flow for separations of the order of 
A. Define the "transverse" structure function across A 
as St (A) = {[ui(Ut) - u 2 (Ut)} 2 } and the "longitudinal" 
structure function as S L ( A) = {[u^Ut+Ut^-u^Ut)} 2 ) 
where t& — A/U. If the flow were isotropic we would ex- 
pect § 



S T (A) = Si (A) + 



A dS L (A) 
2 dA 



(3) 



In the isotropic state both components scale with the 
same exponent, St,l(A) cx A^ 2 , and their ratio is com- 
puted from (§) to be l + C 2 /2 « 1.35 where C 2 ~ 0.69 (see 
below). The experimental ratio was found to be 1.86, in- 
dicating some 40% anisotropy at this scale Q ; we expect 
even more anisotropy on larger scales. 

To obtain a theoretical form of the structure function 
tensor we first select a natural coordinate system. An ob- 
vious choice is the mean- wind direction n along the a — 3 
axis. A second axis is given by the separation vector A 
between the two probes. It turns out that for this par- 
ticular geometrical configuration such a two-dimensional 
coordinate system is sufficient to describe all the non-zero 
tensor components. This is to be expected since all mea- 
surements of velocity and all separations are contained 
in the plane (n, A). We can envision a more general 
experimental setup which measures velocity components 
perpendicular to the ground or in which the probes are 
separated in the vertical direction. In such a situation, 
we would need to take into consideration tensor com- 
ponents arising from the existence of a non-trivial third 
direction orthogonal to the (n, A) plane. We expect that 
such experiments will in the future allow us to study the 
antisymmetric components of the tensor (see below). 

To continue, we need to write down the tensor form for 
the general second order structure function (defined by 
Eq. (^]) for n = 2) in terms of irreducible representations 
of the SO(3) rotation group. This tensor can be written 
in terms of the representations of the direct product of 
two three-dimensional Euclidean vector spaces (for the 
indices a, (3) and the space of continuous functions on 



the unit sphere (for the direction of R) g . The latter is 
spanned by the spherical harmonics Y/. m , and the repre- 
sentations of the product space are indexed by j, denot- 
ing a 2j + l dimensional irreducible representation. Every 
such representation is associated with a scalar function 
Cj(R) which is expected to scale with a universal expo- 
(7) 

nent Q ! the latter is an increasing function of j and 
Q = (2- Previous theoretical considerations |)| led to 
the estimates ( { 2 1} « 1, cf } ~ 4/3. We are interested in 
relatively modest anisotropies, and so focus on the lowest 
order correction to the isotropic (j — 0) contribution. In 
other words, we write 



S Q/5 (R) = Sf£ (R) 



Sff 2 (R) 



(4) 



We do not have a j = 1 term since the possible contribu- 
tions to it vanish either because of parity considerations 
(the structure function itself is even in R) or by the in- 
compressibility constraint. Now the most general form of 
the tensor can be written down by inspection. The case 
j = is well known, and we write it as 



SK J (R) = c (fl) 



(2 + C 2 )<5 Q/3 - C: 



R a RP 
~R 2 ~ 



(5) 



where cq(R) = cqR^ 2 , and cq is a non- universal numerical 
constant that needs to be fitted to the data. The j = 2 
component can be written as 



S^ 2 (R) = c 2 (R) 
R a RP(n ■ R) 2 



a5 *e +b ^ +d ^:*) 2 (6) 



n 2 



R 2 



+ fn « n P + g( "„f ) (R a nP + R?n a ) 
2R Z 



where c 2 (i?) = R 1 * 2 ■ Here a, b, d, e, / and g are un- 
known coefficients which are independent of R. This form 
can be further reduced by imposing the conditions of in- 
compressibility and orthogonality with the j = part 
of the tensor. This leaves us with only two independent 
coefficients of the form 



S°f 2 (R) =aR C2 



(2) 



2)^-Cf(d 2) +6) 



x «$■ 



^(n.fi) 2 +2C ( 2)(c ( 2) 



R 2 



2) 



R a R f3 (n-R) 2 

^4 



ad 2) ] 2 
d 2) (d 2) 



3d 2) 

2), 



6)n a n 



(R a n + R p n a ){n ■ R 



R 2 



(7) 



bR& 



(2) 



3)(C£° +2)5 al3 {n-R) 2 



R a R? 



(Cr+3)(Cr +2)« a n 
R a RP{n-R) 2 



42) 



i? 4 



ac 2 



R 2 

(2d 2) + i)(d 2) -2) 

(2) i 2 - 4)(i? < V + R p n a ){n ■ R) 



Finally, we note that in the present experimental set-up 
only the component of the velocity in the direction of n 
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TABLE II. The scaling exponents and the three coefficients 
in units of (m/sec) as determined from the nonlinear fit of 
Eqs. (9) to the data. 



is measured. In the coordinate system chosen above we 
can read from (Q) the relevant component as 



S 33 (R,t 



R 

c ° ' A 



(8) 



2 + C2 - C2 cos 2 



+a 



+b 



(2) 



d 2) (3Cf ) +2)cos 2 



+ 2Cf ) (Cf -2)cos 4 



(d 2) + 2)(d 2) + 3) - Cf } (3d 2) + 4) cos 2 



+ (2^ + 1)^ 



2) cos 4 9 



Here 9 is the angle between R and n, and we normal- 
ized R by A making all our coefficients dimensional, with 
units of (m/sec) 2 . 

To fit our experimental results we converted, using 
Taylor's hypothesis, the structure functions computed 
from time differences for a single probe, and cross dif- 
ferences between the two probes, to components of the 
form (|^) with 9 — and with variable 9, respectively. In 
other words, 

S 33 (R, 6 = 0) = {[u x {Ut + Ut R ) - Ul {Ut)} 2 ) , (9) 
where t R = R/U, and 



S 33 (R, 9) = {[ui{Ut + Ut k ) - u 2 {Ut)] 2 ) 



(10) 



Here 



= arctan(A/[/£^), — R/U, and R 



2 + (Utji) 2 . We then found the three unknown coef- 

(2) 

ficients Co, a and b, and the two exponents C2 and Q 
by a nonlinear fit to the form mj. The fitted numbers 
are presented in Table 2. 

One surprise of the analysis is that the 9 = (purely 
longitudinal) structure function is hardly affected by the 
anisotropy. The reason is the very close numerical abso- 
lute values of the coefficients a and b. In the case 9 = 
the two tensor forms multiplied by a and b coincide, and 
the j = 2 contribution becomes very small. We believe 
that this is the main reason that the effect of anisotropy 
has been disregarded by and large in atmospheric exper- 
iments. The component S 33 (R, 9 = 0) can be reasonably 
fitted by a power law with single exponent. In fact, the 
value of (2 = 0.69 quoted above can be obtained from 
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FIG. 1. The structure functions S 33 for 
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in panels 

(a) and (b), and for non-zero 8 in panels (c) and (d). The 
dots are for experimental data and the line is the analytic 
fit. Panels (a) and (c) present fits to the j = component 
only, and panels (b) and (d) to components j = and j = 2 
together. 
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FIG. 2. The determination of the exponent Q from a 
least-square fit of S 3i (R, 0) to its analytic form. The numer- 
ical value of the best fit, £ = 1.36 ± 0.1 is in close agreement 
with the theoretical expectation of 4/3 (without intermittency 
corrections) . 



merits) are just the smallest exponents in the hierarchy 
Q that characterizes higher order irreducible represen- 
tations indexed by j. The study of these exponents is in 
its infancy, and considerable experimental and theoretical 
effort is needed to reach firm conclusions regarding their 
universality and numerical values. We expect the expo- 
nents to be a non-decreasing function of j, explaining 
why the highest values of j are being peeled off quickly 
when R decreases. Nevertheless, the lower order values 
of can be measured and computed. It is the program 
of the present authors to proceed in this direction. 
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such a simple fit, and as long as one measures only this 
component, there is no reason for a more sophisticated 
analysis. In fact, £2 is so well estimated that we can safely 
take it as fixed when we perform the simultaneous fit of 
all the other unknowns. 

On the other hand, the finite 9 components are 
strongly affected by the anisotropy, and the inclusion 
of the second exponent Q is essential for a good fit. 
In Fig. 1 we exhibit various quantities as a function of 
R/A; the first panel shows S 33 (R 7 9 = 0) with its best 
fit to the j — contributions, and in panel (b) with its 
best fit to the sum of j = and j = 2. The differences, 
though small as noted previously, are enough to produce 
an improved agreement. Similarly in panels (c) and (d) 
we show S 33 (R, 9) with its best fit to the j = and the 
sum of j = and j = 2, respectively. It is clear that 
panel (c) does not give an adequate fit in the "inertial 
range" , as estimated from panel (a) to range up to, say, 
15A. If we tried to read a log-log slope from the data in 
panel (c) we would get an apparent scaling exponent £2 
considerably smaller than 0.69. The excellent fits in pan- 
els (b) and (d) are a good support to the present mode 
of analysis. 

(2) 

To our knowledge, this determination of Q is the first 
instance of finding an exponent describing the degree of 
anisotropy. The close agreement with the theoretical ex- 
pectation of 4/3 (e.g., Ref. is a strong indication 
that this exponent is universal. In Fig. 2 we present \ 2 
(the sum of the squares of the differences between the 
experimental data and the fitted values) as a function of 
Q • The optimal value of this exponent and the unccr- 

(2) 

tainty determined from this plot is Q ~ 1-36 ± 0.1. 

It should be understood that the exponent Q (and 
also £2 that is unavailable from the present measure- 
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